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ABSTRACT 

Relativistic magnetized jets are key elements in Active Galactic Nuclei and in other as- 
trophysical environments. Their structure and evolution involves a complex nonlinear 
physics that can be approached by numerical studies only. Still, owing to a number of 
challenging computational aspects, only a few numerical investigations have been un- 
dertaken so far. In this paper, we present high-resolution three dimensional numerical 
simulations of relativistic magnetized jets carrying an initially toroidal magnetic field. 
The presence of a substantial toroidal component of the field is nowadays commonly 
invoked and held responsible for the process of jet acceleration and coUimation. We 
find that the typical nose cone structures, commonly observed in axisymmetric two- 
dimensional simulations, arc not produced in the 3D case. Rather, the toroidal field 
gives rise to strong current driven kink instabilities leading to jet wiggling. However, 
it appears to be able to maintain an highly relativistic spine along its full length. By 
comparing low and high resolution simulations, we emphasize the impact of resolution 
on the jet dynamical properties. 
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1 INTRODUCTION 

Understanding the dynamics of relativistic jets is an es- 
sential step for the comprehension and interpretation of 
many phenomenological properties of Active Galactic Nu- 
clei. Magnetic fields appear to be an essential ingredient in 
their structure: observations of highly polarized non-thermal 
synchrotron radiation evi nce the presence of p artially or- 
dered magnetic fields, e . g.jGabuzda et al (|2004l ) ; iLaing et all 
(|2008l '): [Ga"buzda et all (|2008l ). In particular, it is Hkely that 
a toroidal component prevail on large scales as the poloidal 
component decays as Bp oc 1/-Rj while the toroidal oc 
1/Rj, where 1/Rj is the jet radius that expands from the 
inner galactic core to the radio jet lobes. From the theoret- 
ical point of view toroidal magnetic fields are the necessary 
element of the widely accepted magnetically driven mecha- 
nisms for jet acceleration and coUimation by fields with foot- 
loints attached to a spinning b l ack hole or accretion disc 



Shafranov criterion (|Batemanlll978l ) 



Bj, 



Br, 



27r7?, 



where z is the jet length. Instead astrophysical jets appear 
to be quite stable and this becomes an unsolved issue, al- 
though several proposals have been made, related to the 
for mation of cocoon s by the bow shock at t he head of the 
jet (lBegeImanlll998|'). jMassaglia et al.lll996f) or the jet ex- 



pomts attacnea to a spmnmg p l acJt nole or accretion disc 
(lLovelacdll976l: iBlandfordll 19761 : iBlandford fc Pavne|[l983 : 
iNaravan et al.l^2007^ . 



On the other hand cylindrical magnetohydrodynamic 
(MHD) configurations are subject to many unstable modes, 
as reflection modes, Kelvin-Helmholtz modes, current carry- 
ing modes, etc. In particular configurations in which toroidal 
magnetic fields dominate are known to be violently unsta- 
ble to the m = 1 kink instability according to the Kruskal- 
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pansion (|Mo11 et al.ll2008l : iMcKinnev fc Blandfordllioogj '). or 
wave coupling with MHD modes producing particle acceler- 
ation (jBenford et al.lll980l ). 

On the other hand, one may wonder whether the main 
jet constituent is electromagnetic energy in the form of 
Poynting-dominate d beams, rather than mass. However, 
ISikora et al.l (|2005l ) argued that, even though jets could be 
Poynting-dominated at the origin, observational data sug- 
gest that they become kinetically dominated beyond about 
1,000 gravitational radii from the central acceleration region, 
powered by a super-massive black-hole of ~ 10*~^° M©. 

The mechanism by which the toroidal field responsi- 
ble for the jet acceleration can be dissipated may be re- 
lated to the above ment ioned instabiltics. In this respect 
Giannios fc SpruitI (|2006l ) (see also ,Eichlcrill993l : iBegelmanI 
1998h analyzed the role of the kink instability in magnet- 



ically driven jets and concluded that the Poynting flux 
can be rapidly dissipated and the jet becomes kinetically- 
dominated, again at about 1,000 gravitational radii. 
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Recently, iMcKinnev fc BlandfordI (|2009l ) showed that 
relativistic jets produced in proximity of rapidly spinning 
black-holes with dipolar magnetic fields expand rapidly and 
can survive up to large radii without showing appreciable 
disruption due to the dominant m = 1 kink mode. Their 
simulations were carried out using large-scale global gen- 
eral relativistic simulations in three-dimensional spherical 
coordinates adopting a modest resolution to capture both 
the disk and the jet dynamics. Three-dimensional simu- 
lations of the jet formation and evolution in the general 
relativistic MHD dom ain has been performed as well by 
iNishikawa et al.1 (j2005h in the case of a Schwarzschild black 
hole accretion disk system. 

It is therefore crucial to study in a more systematic fash- 
ion the fate of relativistic magnetized jets in this transition 
region, examine the model parameter space to find whether 
the above requirements and predictions can be met. Rel- 
ativistic jets carrying magnetic fields can be described by 
the equations of relativistic MHD (RMHD) and only in the 
last decade considerable efforts have been spent towards ob- 
taining more robust and accurate numerical schemes. This 
opened the door to noticeable advancements in the field 
of numerical simulations. A first study of the properties 
of two-dimens ional relativ is tic m agnetized jets has been 
performed by iKomissarovl (|l999l ) and subsequent l y the y 
have been thoroughly analyzed by Leismann et al.l (l2005h: 
iMignone eraP (j2005l '): iKeppens et al.1 (|2008l ) . iMizuno et al.1 
have, instead, studied the instability of a 3D RMHD 
cylindrical flow. 

In this paper we present the flrst high-resolution (to 
our knowledge) numerical simulations of the propagation in 
three dimensions of a (special) relativistic magnetized jet in- 
jected in a uniform un-magnetized external medium. Com- 
putations are carried out at both high (640 x 1600 x 640 
ceU, with 20 pts/beam radius) and low (320 x 800 x 320, 
with 10 pts/beam) resolution. A comparison between the 
two reveals the importance of high resolution, especially 
for this class of astrophysical phenomena. We consider jets 
initially carrying a purely toroidal magnetic field; as pre- 
viously said, the presence of this component is intimately 
connected to the jet formation mechanism and it is likely to 
affect the jet dynamics through the onset of current driven 
(CD) instabilities. In this respect, it is of great interest to 
analyse their growth and evolutio n. Previous non-rela t ivistic 
three dimensional simulations of iNakamura fc Meieil ()2004l l 
revealed that strongly magnetized supersonic jets develop 
CD induced distortions and that the kink (m = 1) mode 
;rows faster than the other modes. More recently, IMoU et al.l 
20081 ) have studied the generation and propagation of a 
nonrelativistic jet with a self-consistent toroidal field com- 
ponent, showing and analysing the development of these in- 
stabilities. In the present work we investigate the relativistic 
case, starting with the limiting case of a purely toroidal field 
which, at least in the Newtonian case, is known to be the 
most unstable configuration. Our study extends previous 3D 
simu lations of the pro pagation of unmagnetized relativistic 
jets ([Rossi et al.ll200^ ) and, in the following, we will often 
refer to these results for comparison. The effects of different 
field configurations and physical parameters will be exam- 
ined in a forthcoming paper. 

In the next section we describe our numerical setup, in 



^we outline the simulation results and in ^we summarize 
our conclusions. 



2 NUMERICAL SETUP 

Numerical simulations are carried out by solving the equa- 
tions of spec ial r elativistic MHD (RMHD) in conservation 
form fiAnileiilS)90l) : 



9m (p^*" 



= 0, 



(1) 



0, 



where p is the rest mass density, = 7(1, u) and 6'' — 
(fo°,i?/7 -I- b"v) are the four-velocity and covariant mag- 
netic field written in terms of the three velocity v and lab- 
oratory magnetic field B and b° = ■ B. A flat metric 
^i^v _ (jiag(|_x, 1, 1, 1) is considered. In the previous equa- 
tions w — vjg_ + and p = pg + }? jl express the total 
enthalpy and total pressure in terms of their thermal (wg 
and pg) and magnetic contributions (\? = b'^b^), respec- 
tively. An additional equation, describing the advection of 
a passive scalar (or tracer) f{x,y,z,t), is included to dis- 
criminate between jet material (where / = 1) and the envi- 
ronment (where / = 0). We assume a single-specie relativis- 
tic perfect fluid (the Synge gas) des cribed by the approxi - 
mat ed equation of state proposed bv lMignone et al.1 (|2005l ) 
and iMignone fc McKinnevI (120071 ) . 

The computational domain is initially filled with an un- 
magnetized medium at rest, with uniform density pa and 
gas pressure pa- The jet fiows parallel to the y direction 
through a cylindrical nozzle specified by r = \/x'^ + < 1 
(in units of the jet radius) at the lower boundary. At the 
nozzle, the beam has constant density pj and longitudinal 
velocity component (vy) given in terms of the Lorentz factor 
7j, whereas the magnetic field carries a purely toroidal (i.e. 
azimuthal) component: 



Bmr/a for r < a , 
B4,{r) = { Bma/r for r>a. 



(2) 







for r > 1 . 



where a = 1/2 is the magnetization radius. The pressure 
profile follows from the numerical integration of the radial 
momentum balance equation across the beam: 



d , 



(3) 



where a small differential rotation is assumed by setting 
= ayBtp/ Bm with a — 0.2. The pressure profile p{r) 
is continuous across the jet boundary, i.e., pa ~ Pj(l)- The 
integration constant and the field strength B^ are speci- 
fied through the sonic Mach number Ala = Vj/ca{a) and 
Alfvenic Mach number Ma = Vj/vA{<i) {cs and va are the 
sound and Alfven velocities) at the magnetization radius. 
Thus the beam values are entirely specified in terms of the 
density contrast 77 — pj/pa, Lorentz factor 7^, sonic and 
A lfvenic Mach numbers. Motivated by the results obtained 
in lRossi et all (|2008l ). we choose rj = 10"*, ■yj = 10, Ms = 3 
and Ma = 1.67, so that the ratio of the average gas and 
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Figure 2. Volume renderings of tiie tliermal pressure distribu- 
tions at t = 710 for the high resolution 3D run. 



magnetic pressures is close to one (equ ipartition). Velocit y 
perturbations are super-imposed as in I Rossi etaL I (|2008D . 
Reflective boundary conditions hold at y = outside the 
inlet region and zero-gradient is assumed through the re- 
maining boundaries. 

Simulations are carried out with the PLUT O code for 
astrophysical gasdynamics (jMignone et al.1 1200/1 '). The cho- 
sen configuration employs a second order Godunov type 
scheme with the rece ntly developed HLLD Riemann solver 
(|Mignone et al.l [20091 ) and the mixed hyper bolic / parabolic 
divergence cleaning by iDedner et al.l (|2002l ') to control the 
solenoidal constraint. Our computational box is defined by 
X G [-L/2,L/2], y G [0, Ly] and z G [-i/2,L/2], (i = 56, 
Ly = 80) with a uniform resolution in |a;| ,\z\ < 10 and a geo- 
metrically stretched grid otherwise. Both low (10 points per 
beam radius) and high (20 zones) resolution simulations are 
performed, yielding an overall resolution of 320 x 800 x 320 
and 640 x 1600 x 640 computational zones, respectively. 



3 RESULTS 

In the left panel of Fig. [T] we show the volume renderings of 
the tracer distribution when the jet has reached the length of 
~ 70 beam radii. For the sake of comparison we also show, in 
the right panel, the rendering of an equivalent 2-D axisym- 
metric simulation with the same parameters. The picture 
(in the right panel) clearly indicates the presence of a nose 
cone structure typical o f two-dime nsional high Poyriting fiux 
jets (IClarke et al.lll986l : iLind et al. .1989; Komis sar ov 1999i ; 
iLeismann et al. 20051 ). originating from the amplification at 



the terminal shock of the toroidal field component which 
then confines the jet material and prevents it to freely flow 
in the cocoon. In three dimensions, however, this structure 
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Figure 3. Barycenter Displacements from the longitudinal axis f 
of the Lorcntz factor (top) and jet laboratory density (bottom) for 
the high resolution (solid lines), low resolution (dashed lines) and 
hydro (dotted lines) cases. The corresponding integration times 
are t = 710, t = 840 and t = 590, respectively. 



becomes unstable and does not form leading to a very dif- 
ferent asymmetric morphology. In following the jet propa- 
gation, we observe that its trajectory becomes progressively 
more curved, moving away from the longitudinal y axis. This 
effect becomes more pronounced at the jet head and can be 
attributed to the presence of current driven kink instabili- 
ties. Moreover, the backflowing material forms a very asym- 
metric cocoon as a result of the changes in direction of the 
jet head. In Fig. [21 representing the pressure distribution, 
it is possible to observe two regions of enhanced pressure 
in the bow shock. Comparing the tracer and pressure dis- 
tributions, we can identify these areas as the points where 
the jet is more strongly deflected. The asymmetric backflow 
appears to be quite pronounced, reaching relatively high ve- 
locities up to 0.9c. 

The displacement of the beam from the longitudinal 
axis plays a fundamental role in the jet morphology and it 
is evaluated, at each position along the jet, by the quantity 
f(y) — \/x^ + z^, where 



x{y) = 



J xQdx dz 



4y) = 



J zQdx dz 



(4) 



J Qdx dz ' J Qdx dz 

where Q can be any flow quantity. Choosing Q — 7X, where 
X = 1 for 7 > 1.5 and x ~ ^ otherwise, gives a measure of the 
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Figure 1. Volume renderings of the passive scalar distributions for the high resolution 3D run (left panel) and 2D axisymmetric case 
(right) at t = 710. 



distance of the Lorentz factor centroid from the axis. Like- 
wise, we can quantify the jet density barycenter by choosing 
Q = p7/. In Fig. [3] we plot both the Lorentz factor and 
density centroids as functions of y. Sohd, dashed and dotted 
hnes refer, respectively, to the high resolution, low resolu- 
tion and non magnetized case with the same parameters 
(Lorentz factor, density ratio and Mach numbers). Simula- 
tion times have been chosen so that the jets have reached 
approximately the same distance, ~ 70 radii. Both the low 
and high resolution cases show moderate displacements up 
to ~ 40 jet radii followed by a strong increase soon after. 
The effect is more pronounced in the high resolution case, 
eventually reaching maximum values in excess of 10. On the 
contrary, no significant deviations are seen in the purely hy- 
dro case and the beam propagates very close to the initial 
axis. 

Fig.|3]represents the distribution of the jet mass fraction 
moving at a certain value of 7/3 at the end of the simula- 
tion when the jet has traveled approximately 70 beam radii. 
Both low (top panel) and high (bottom panel) resolution 
cases are shown. We use the four-velocity (instead of the 
Lorentz factor) to avoid compression of the scale close to 
7 = 1, i.e. at low velocities. The legend on the right gives 
the corresponding value of 7/3 for each color. As an illustra- 
tive example, one can see that at y = 15 the mass fractions 
moving with 7/3 > 6, 7/3 £ [5,6],7/3 £ [4,5],7/3 e [3,4] and 
7/3 < 3 are, respectively, 10%, 10%, 20%, 25% and 35%. 
Focusing on the high velocity part of the distribution, we 
can observe the presence of material moving at 7 ~ 10 all 
along the jet (red spots close to the y-axis) . The behavior of 
the low resolution case, shown in the upper panel, is quite 
different. Jet material moves at 7 ~ 10 only up to j/ ~ 20, 
with a decrease of the maximum Lorentz factor to a value 
of about 5. 



Mass Fraction, t = 840,0 7/? 
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Figure 4. Jet mass fraction for the low resolution (top panel) 
and high resolution (bottom panel) cases. 

The toroidal field seems to protect the jet core from any 
interaction with the surrounding, thus preventing momen- 
tum transfer to the external medium. This behavior is bet- 
ter understood by inspecting the top panel in Fig. [5] where 
we show the jet cross-section (in the high resolution case) 
at different positions along the longitudinal axis. For com- 
parison, a similar plot for the unmagnetized case, shown 
in the bottom panel, reveals that in this case the jet sur- 
face is deformed by non-axisymmetric modes of high order 
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Figure 5. Two dimensional slices at constant y planes for the high resolution magnetized case (top panel) and purely hydrodynamical 
case (bottom, low resolution). The solid line marks the longitudinal axis (x, 2 = 0) while the line of sight is such that the jet propagates 
from the furthest boundary to the closer one. 



(m > 2) promoting mixing between jet material and external 
medium. On the contrary, the high resolution magnetized 
case shows the presence of low order modes (m ^ 2) only, 
producing oval deformations that cannot induce any mixing. 
The tension of the toroidal field, in fact, acts as a strong sta- 
bilizing factor for high o rder modes. Our re sults are thus in 
agreement with those of iRoseri et al.1 (| 19991 ) where the pres- 
ence of a toroidal field was shown to noticeably lessen the 
growth and amplitude of small wavelength modes. 

From Fig.|3]one notices also the periodic decrease of the 
maximum Lorentz factor associated with jet pinching and 
the formation of periodic shocks along the beam, also visi- 
ble in the left panel of Fig. [51 where it is possible to follow 
the appearance of periodic compressions all along the jet. 
The formation of these internal shocks has been observed 
in all jet simulations and is in general attributed to the in- 
teraction of the cocoon with the jet proper and to Kelvin- 
Helmholtz instabilities. Three-dimensional effects seemed to 



decrease the importa nce of these intern al shocks. Indeed, in 
a similar hydro case (|Rossi et al.ll200^ ) at the same resolu- 
tion, we observed their presence only in the first half of the 
jet. Here, the decreased interaction of the jet core with the 
ambient favors their formation and, besides, the additional 
pinching effect of the toroidal field increases the compres- 
sion factor. At the same time, the presence of the toroidal 
field has also the effect of reducing the distance between the 
shocks when compared to the hydro case. We have seen that, 
in the case that we are examining, the jet core is protected 
by the toroidal field and maintains its initial Lorentz factor 
all along its length. However, it is interesting to ask whether 
an effect of mixing and entrainment can actually take place 
on the more external jet layers. From Fig. 3] we can observe 
that the mass fraction of the material with 7/? < 2 varies 
between 20% in the region closer to the jet origin, to about 
30% at outer positions, while in the hy dro case th i s frac - 
tion reached values of more than 50% (|Rossi et al.ll2008l ). 



6 Mignone et al. 



showing a smaller entrainment efficiency in this case. Cor- 
respondingly, a larger fraction of the jet material moves at 
higher values of 7/3. 

These results show that, in spite of the large deflec- 
tion induced by the presence of strong toroidal magnetic 
field, non-axisymmetric m = 1 screw modes do not lead 
to appreciable je t disruption. These findings agree with 
those reported in iMcKinnev fc BlandfordI (|2009l ) where the 
m = 1 modes are effective only in modeling the jet substruc- 
ture without affecting the large-scale dynamics. Although 
our simulation setup is noticeably different, this qualitative 
agreement supports the general idea that some other inter- 
vening mechanism must act so as to suppress t he instability 
growth, (|Mo11 et al.ll2008l : iNaravan et al.ll2009l '). 

The simulation we are discussing has been performed 
at high resolution employing 20 grid points over the beam 
radius. Numerical resolution effects are known to be impor- 
tant in jet simulations and, for this reason, we have com- 
pared the high resolution computation with a lower resolu- 
tion one, with only 10 points on the jet radius. In the low 
resolution case, as it is shown in the upper panel of Fig. [l] 
numerical diffusion effects tend to increase the deposition of 
momentum from the jet core to the ambient and therefore 
the jet is not able to maintain its initial Lorentz factor up to 
the head. In addition, we have also observed that periodic 
shocks along the jet appear to be substantially weaker and 
they are not present all along the jet as in the high resolu- 
tion case. Besides, as discussed above, the jet displacement 
results larger at high resolution. This comparison shows that 
the differences due to grid resolution are substantial and 
demonstrates the importance of performing high resolution 
simulations for reproducing more accurately the physical be- 
havior. A question that may occur is obviously whether the 
resolution employed in this work is sufficient, since our re- 
sults show that we have not yet reached convergence in the 
basic jet properties (a complete convergence for this prob- 
lem cannot be reached). We have shown a tendency towards 
larger displacements of the jet off the axis as the resolution 
is increased, although this effect is, in any case, localized 
at the jet's head. Still, we observe a reduction of the mo- 
mentum transfer from the inner jet to the ambient as the 
numerical resolution increases. We expect that these well de- 
fined trends in the physical characters of the solution should 
be maintained by a further increase in resolution, that may 
be feasible in future studies on this subject. The effects de- 
scribed above may therefore change quantitatively, but the 
qualitative features that we outlined should remain similar 
to what we obtain. 

The results pertaining to the hydrodynamical cases dis- 
cussed in I Rossi et al.l (|2008l ) indicated the density ratio 77 as 
the main parameter governing the interaction between the 
beam and the external medium. In order to investigate the 
role of 77 for the present magnetized model, we performed 
a low resolution computation with rj — 100. The results 
show that, up to the same length, the jet propagates al- 
most undisturbed, without any sign of kink instability nor 
of mixing with the ambient medium. We think that the most 
likely explanation for the difference between the two cases at 
different density ratios (this parameter has no direct effect 
on the kink instability) should be related to the different 
advance velocity of the jet head, that is much larger for 
?7 = 100 than for — 10, 000. In fact, the unstable perturba- 



tions are advected by the jet as they grow. In the case with 
rj — 10, 000 the jet head advances at a velocity that is much 
lower than the jet velocity and, therefore, the unstable per- 
turbations accumulate at the jet head and have enough time 
for their growth and can give origin to the jet deformation. 
In the lower ratio case, instead, the jet head advances at a 
much higher velocity and the kink instability does not have 
enough time to grow. In this case its effect should be visible 
only on a larger scale. 



4 CONCLUSIONS 

We have presented the first results of a three-dimensional 
high-resolution numerical study of the propagation of rel- 
ativistic magnetized jets injected into a uniform medium 
using the PLUTO code. At the highest resolution, the com- 
putations employed 20 points on the jet radius for an overall 
resolution of 640 x 1600 x 640 zones covering a Cartesian do- 
main of 56 X 80 X 56, in units of the beam radius. 

From the variety of possible magnetic field configura- 
tions, we decided to start our study by considering jets ini- 
tially carrying a purely toroidal magnetic field. Indeed, the 
most accepted jet formation mechanisms predict that this 
field component has to be present and this strengthen the 
necessity of understanding its key role on the jet dynamics. 
The presence of a toroidal component of the field is known 
to drive current-driven kink instabilities which in our case 
are responsible for jet wiggling and beam d efiection off the 
main longitudinal axis (|Naravan et al.ll2009l ') . The deviations 
tend to become more prominent towards the final parts of 
the jet, where we measured off-axis displacements in excess 
of 10 radii over a length of 70 radii. The wandering of the jet 
head, induced by kink instability effects, may create multiple 
sites where the jet impacts on the external medium form- 
ing strong shocks. This behavior may originate the multiple 
hot spots that are observed in several rad iogalaxies (see e.g. 
iLeahv et aLlll997l : [Hardcastle et al.ll2007^ . The asymmetry 
of the backflow is another feature that replicates the obser- 
vational appearance of several objects. 

An additional effect of the toroidal field component 
seems that of shielding the inner jet core from any inter- 
action with the surroundings and therefore from a loss of 
momentum. Furthermore, we have observed that the jet re- 
mains highly relativistic all along its length. Quasi periodic 
shocks, formed by the combination of several effects like the 
interaction with the cocoon and the toroidal field pinching, 
can be observed all along the jet extension. A spine plus 
sheath layer structure with shocks supports the idea that jets 
at sub-kpc scale can originate correlated variability at radio 
frequencies and in X and gamma rays (see for instance recent 
observations by t he AGILE and FERM I missi ons and the 
TeV g round arravs lPonnarumma et al] (120091: Acciari et al.l 
Hooi)). The very high frequency emission would originate 
in the slower sheath layer by inverse Compton boosting of 
the optical and X photons produced in the relativistic spine. 
The X ray emission would be synchrotron emission not self- 
absorbed and would have variability strictly connected to 
the gamma ray emission. The radio emission would come 
from s ynchrotron emission in an o uter region of the jet 
sheath (iTavecchio fc GhisellinJ|2008r ). 

Recently Moll et al.l ( 200S ) have performed 3D Simula- 
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tions of a Newtonian magnetized jet with a toroidal compo- 
nent and have shown the development of kink instabilities 
that produce helical jet deformations. In their setup the jet 
is generated by maintaining a rotation profile in the bound- 
ary zones and a spherical grid, with a non uniform spacing in 
the radial direction, is used to describe the jet propagation. 

The comparison between low and high resolutions cases 
shows significant differences in the formation of shocks along 
the beam. In the high resolution simulations, shocks appear 
to be stronger, closer and in a larger number. In addition, 
the jet is able to keep a highly relativistic core for a longer 
distance. 

The results presented here represent only a first step in 
the analysis of the behavior of relativistic magnetized jets. 
We have shown that three dimensional and high resolution 
effects are important factors in this investigation. In a forth- 
coming paper we intend to explore the infiuence of different 
(e.g. poloidal or helical) field configurations. 
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